" UP: Have Observations Revealed a Variable Upper End of the Initial Mass Function ? " 
ASP Conference Series, Vol. 4XX 

Marie Treyer, Ted Wyder, Don Neill, Mark Seibert, Janice Lee, eds. 
©20XX Astronomical Society of the Pacific 



Measuring the Upper End of the Initial Mass Function with 
Supernovae 



James D. Neill 1 

1 California Institute of Technology, 1200 East California Blvd., Pasadena, CA 
91125 

Abstract. Supernovae arise from progenitor stars occupying the upper end of the 
initial mass function. Their extreme brightness allows individual massive stars to be 
detected at cosmic distances, lending supernovae great potential as tracers of the upper 
end of the IMF and its evolution. Exploiting this potential requires progress in many 
areas of supernova science. These include understanding the progenitor masses that 
produce various types of supernovae and accurately characterizing the supernova out- 
burst and the environment in which it was produced. I present some preliminary work 
identifying the environmental conditions that produce the most luminous supernovae, 
believed to arise from stars with masses greater than 100 M . I illustrate that the pres- 
ence of these extreme supernovae in small star-forming dwarfs can be used to test our 
understanding of the upper end of the IMF. 



1. Motivation 

Supernovae (SNe) are the spectacular end result of massive star evoluti on. Current 



estim ates of the lower mass limit for SN progenitors are near 8 M Q (e.g. JSmartt et al. 



120091) . making them extremely relevant to studies of the upper end of the initial mass 



funct ion (IMF). The fact that SNe are now being detected beyond z = 2 (ICooke et al. 
2009) adds the potential to use SNe as tracers of IMF evolution. Their usefulness 
depends on how well their progenitors can be understood and how well the observations 
of SNe can constrain the distribution of progenitor masses. In addition, we must be able 
to characterize the host galaxies of SNe and connect them to our general understanding 
of star formation in galaxies at various stages of their evolution. 



2. SN Wish List 

Let us take an idealized view of SNe and list what they have to offer for studies of 
the upper-IMF. For a given galaxy with a given star formation history (SFH) or current 
star formation rate (SFR), we would like to know: what is the mass of the largest star 
that can be produced? This is a fiducia l point in the IMF th at calibrates the relationship 
betwee the IMF and the host SFR (e.g. JWeidner et al1l2010l) . The fact that such massive 
stars are bright is offset by their very short lives, making them difficult to observe. The 
SN explosion marks these massive stars in events detectable across cosmic distances. 
Provided we can decode the progenitor mass and measure the host galaxy SFH/SFR, 
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we can directly calibrate this important fiducial point in a system's IMF, and trace its 
evolution potentially out beyond z — 2. 

If we know the masses of not just the most massive SN progenitors, but a whole 
range of progenitors down to the 8 M Q limit, we can measure the distribution of massive 
stars from 8 to over 100 M Q as a function of host SFR/SFH and mass and, while we are 
at it, metallicity. Assuming we have stellar evolution models for these massive stars that 
will provide us with their lifetimes, we can then look into the evolution of clusters with 
greater accuracy, revealing how energy and processed metals are re-introduced into the 
birth cloud of our idealized SN progenitors. 

Now that we have imagined this ideal universe in which understanding the upper 
end of the IMF rests only on the collection of enough SNe and the corollary host data 
to sample the distributions well enough, let us now look at the reality of SN progenitors 
and how close we can get to this idealized situation. 



3. Reality 



Type 


Mechanism 


Character 


Mass (M Q ) 


Progenitor 


la 


Thermonuclear 


mass-transfering binary 


3-8 


CO WD 


H(P,L,n) 


Core-Collapse 


H in absorption (emission) 


7-25 


R/BSG 


lib 


Core-Collapse 


weak H simlar to He 


15-25 


(late) WN 


lb 


Core-Collapse 


all H removed 


25-40 


(early) WN 


Ic 


Core-Collapse 


all H, He removed 


40-80 


WC/WO 


Iln-lum 


Core-Collapse 


narrow H emission 


80-150 


LBV 


Ipec 


pulsational pair? 


all H, He removed 


80 - 150 


W-R? 


I-PP 


pulsational pair 


Fe-group but no H, He 


>100 


pop III? 


Table 1 


SN Types 









I will briefly outline the menagerie of SN types and what we know about their 
progenitors. I have summarized the basic prop e rties o f SNe of various types in Table 1 
after the scheme presented in iGal-Yam et all (120071) . SNe are divided roughly into 
two categories depending on the presence (type II) or absence (type I) of hydrogen, 
either in emission or absorption, in their spectra. When we attribute physical processes 
to the SN explosions it becomes clear that a more rational division is between those 
explosions caused by the collapse of a massive stellar core (core-collapse, or CC SNe) 
and those explosions caused by the thermonuclear burning of a carbon-oxygen white 
dwarf (thermonuclear, or type la SNe). We will not consider th e thermonuc lear SNe, 
since we know they arise in binary systems (for a review, see lLivioll200lh and it is 
unlikely that their progenitors arise from stars more massive than 8 M Q , otherwise they 
would produce CC SNe. 

Within the CC SN types, there appears to be a natural progression of atmospheric 
depletion due to mass loss of varying efficiency. The least depleted are the type II SNe 
that show hydrogen lines. These SNe are further sub-divided by light curve shape (II-P 
for 'plateau' and II-L for 'linear decline') and the presence of narrow emission lines 
(Il-n), indicating the interaction of the SN ejecta with surrounding circum-stellar mate- 
rial. The type lib SNe represent the next level of depletion in which helium lines are as 
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strong or stronger than the hydrogen lines. Type lb SNe occur when the atmospheric de- 
pletion is so complete that only helium lines are seen. The most depleted objects are the 
type Ic SNe, showing neither hydrogen nor helium lines in their spectra. Some of these 
depleted objects have cores so massive that they undergo a pulsational pair-production 
instability which eventually produces a very luminous explosion. These ar e the pair in- 
stabil i ty SNe (PISNe) that have been theoretically predicted for de cades dBarkat et all 
119671 : iBondet all 119841 : iHeger & Woosleyl 120021: IWaldmanl 120081). but only recently 
confirmed observationally (iGal-Yam et alj |2009: Oui mbv et alJl2009h . 

The cause of the mass loss that depletes the hydrogen and helium is some com- 
bination of low surface gravity, stellar winds, and binary interaction. This leads to 
significant uncertainty when estimating progenitor initial masses. We know that metal- 
licity plays a role in determining the efficiency of wind-driven stellar mass loss, but 
there is significant uncertainty in this theory. The ugly reality is that the type of the SN 
derives from some unknown combination of metallicity, binary interaction, and intitial 
stellar mass, clouding our vision of SNe as ideal probes of the upper IMF. 

As of today, the b est constrained SN progenitors are the least depleted ones near 
the lower mass limit. Sma rtt et all ( 20091) report pre-explosion multi-band photometry 
of at least seven SNe II-P and limits on another 13. These result in initial mass estimates 
for the progenitors in the range from ~ 7 - 25M Q . This is a good beginning, but this 
sample represents a small fraction of the population of SN II SNe. There is no way to 
calculate an initial mass distribution from this small sample. 



4. The Dawning Time-Domain Era 



Survey 


Tel. 


FOV 


Filters 


Cadence 


Depth/day 


Coverage 


PTF 


1.2m 


7° 


gR 


l-5day 


~21 mag 


8000 2o /yr 


PanSTARRS 


1.8m 


7° 


grizy 


4day 


22.3(y), 25.1© 


3;r - 30,000 2 ° 


CSS 


0.7m 


8° 


unfit. 


lday 


~20 mag 


1200 2o /dy 


SkyMapper 


1.3m 


6° 


uvgriz 


4hrs-lyr 


-22(g) 


2;r - 26,000 2 ° 


LSST 


6.5m 


10° 


ugrizy3 


3day 


24.5(r) 


3300 2o /dy 


Table 2. 


Current and Future Transient Surveys 







There is no doubt that progress on the upper IMF requires improvements in the 
theory of high-mass stellar evolution. More precise and divers models of SN explo- 
sions will also help in deciphering progenitor characteristics from the spectra of SNe 
in outburst. The engine capable of driving meaningful theoretical progress is a more 
comprehensive and representative data set covering the full diversity of SN explosions. 
We are on the verge of a quantum jump in the quality and the quantity o f observations 
of S Ne in ou tburst with w ide-field transient surveys like PanSTA RRS dHod app et al. 
120041) . PTF dLaw et alJ2009h . CSS dDrake et al.ll2 009), SkyMapper dKeller et alil2007l) 



already online and others such as LSST dlvezic et al.ll2008h coming online in the future. 

These surveys, summarized in Table [2j take advantage of wide-field detectors on 
robotically controlled telescopes to produce large area time-domain surveys that are 
already having an impact on our understanding of SN demographics. Old SN surveys 
used narrow-field detectors and were neccessarily host-targeted, producing biases tied 
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to the luminosity and mass of hosting galaxies. The new areal searches are not only 
removing these biases, but are discovering new types of SNe. 

These surveys have the potential to vastly improve our measurements of SN de- 
mographics and produce accurate distributions of types and rates of each type. With the 
proper followup, these rates can be corellated with host galaxy SFR, stellar mass, and 
average metallicity. They have already shown their value by discovering a new type of 
extremely luminous SN, but it is likely that other rarities will be discovered as well. 



5. Tantalizing Highlight: Extreme SNe in Extreme Hosts 

The first addition to the SN menagerie made by these areal SN surveys already shows 
the power of comprehensive surveys to illuminate regions of parameter space that have 
been in darkness up to now. Because of the host targeted bias, necessitated by narrow 
field searches, the SNe in faint hosts have been passed over. It turns out t hat the most 
luminous SNe on record are typically found in f aint, small dwarf galaxies (Youn g et al. 



20091 : iGal- Yam et alJl2009l : lOuimbv et al.ll2009h . These luminous SNe (LSNe) include 



very luminous explosions showing narrow emission lines that help to power the later 
light curve and are known as SNe Iln-lum. There is another set of LSNe that show no H 
in their spectra, but have unusual light curves and spectra when compared to any other 
SN type. These we call SNe Ipec. 

Recent analysis shows that two LSNe (SN1999as and SN2007bi) show strong 
evidence of having been produced as the result of the pulsational pair instability process 
dGal-Yam et all l2009h . lence we label them Ic-PR Their light curves are consistent 



with models of extreme LSNe that produce > 10 52 ergs of total energy implying the 
radioactive decay of a large amount (> 4M Q ) of 56 Ni. Their spectra s how the results 
of thi s 56 Ni decay by having a high abundance of Fe-group elements dGal-Yam et al. 



I2009h . These signatures require a very massive initial mass for the progenitor star, 
upwards of 150 M w . There is a lso evidence that the progenitor of SN2007bi was a 
single star ( Gal- Yam et al.ll2009h . It is possible that some of the SNe Ipec that have 



similar properties may also be PISNe, but they may not be as ex treme. 

We analyzed t he hosts of seventeen LSNe dNeill et al.l2010h using GALEX (Marti n et al 



2005h and SDSS (lYork et alj|2000h images to derive their UV-optical colors and esti 
mate their stellar masses (M*) and spe cific SFR (sS FR = S FR/M*). We compare them 
to a larger sample of nearby galaxies dWyder et al.ll2007h to illustrate how extreme the 
hosts are and to determine what phase of galaxy evolution might be demarcated by 
LSNe in such small galaxies. While the incompleteness in the comparison sample at 
the faint end prevents a determination of the relative space density of LSNe producers, 
one can see (Figure [B that LSNe appear to prefer small star-bursting dwarfs. 

The LSN hosts have low total SFR, but high sSFR due to their low masses. The 
progenitor of one of the LSNe (SN2007bi) may exceed 150 M©. Any theory of the 
IMF that posits an upper limit to the most massive star produced based on current SF 
must account for these cases, which may be consi dered limiting. To illustrate this, 



consider Figure 4 of iPflamm-Altenburg et al.l (120071) . Here we see that for SFRs near 
10 _1 M G yr -1 , a 150 M Q star is just barely consistent. This is probably good news for 
such theories, because these object are indeed rare. Currently the errors on both the 
host SFR and the progenitor initial mass are large, but as these shrink, the LSNe and 
their hosts could drive a re-calibration of IMF theories that predict the mass of the most 
massive star produced based on current star formation. 
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Figure 1. Gala xy NUV - r versus M r CMD with extreme SNe indicated 
(Neill et al.ll2010h . SN types are indicated by the symbols with Iln-lum SNe indi- 
cated by squares, Ic-PP SNe (PISNe) indicated by triangles, and Ipec SNe indicated 
by diamonds. Right pointing arrows indicate which half -plane the particular LSN 
host is limited to by the UV-optical photometry, upward pointing arrows indicate 
that the NUV - r color must be redder than the symbol, while the double arrows for 
PTF09cnd indicate the quarter-plane that its host is limited to. The contours repre- 
sent the galaxy density per bin (0.5 by 0.5 mag) of the co mparison sample wh ich is 
derived from the GALEX - SDSS cross-match presented in Wy der et al.l 12007). The 
darkest contour indicates a density of 1056 galaxies per bin, while the lightest bin is 
at 132 galaxies per bin. Below the minimum contour galaxies are plotted individu- 
ally as small points. The LSN hosts appear to favor hosts on the blue side of the blue 
cloud, toward the low-luminosity end. 
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6. Conclusions 

While the reality of using SNe to map out the entire upper end of the IMF appears to be 
a long way off, there is cause for hope for the future. Large area surveys will produce 
unbiased samples of SNe, allowing us to test and improve theories of SN progenitors 
and their production. Even though we may never be able to map out the distribution of 
high mass stars with SNe, it does seem likely that LSNe may arise from some of the 
most massive stars known, providing a crucial calibration point in IMF theory. The fact 
that LSNe appear to favor small galaxies with low mass, low intrinsic SFR, but high 
sSFR makes them even more useful to advancing IMF theory. 
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